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ABSTRACT 

Previous  stressed  oxidation  tests  of  C/SiC  composites  at  elevated  temperatures  (350°C  to 
1 500°C)  and  sustained  stresses  (69  MPa  and  172  MPa)  have  led  to  the  development  of  a  finite 
difference  cracked  matrix  model.  The  times  to  failure  in  the  samples  suggest  oxidation  occurred  in 
two  kinetic  regimes  defined  by  the  rate  controlling  mechanisms  (i.e.  diffusion  controlled  and 
reaction  controlled  kinetics).  Microstructural  analysis  revealed  preferential  oxidation  along  as- 
fabricated,  matrix  microcracks  and  also  suggested  two  regimes  of  oxidation  kinetics  dependent  on 
the  oxidation  temperature.  Based  on  experimental  results,  observation,  and  theory,  a  finite 
difference  model  was  developed.  The  model  simulates  the  diffusion  of  oxygen  into  a  matrix  crack 
bridged  by  carbon  fibers.  The  model  facilitates  the  study  of  the  relative  importance  of 
temperature,  the  reaction  rate  constant,  and  the  diffusion  coefficient  on  the  overall  oxidation 
kinetics. 

INTRODUCTION 

Ceramic  matrix  composites,  specifically  continuous  carbon  fibers  in  a  silicon  carbide 
matrix,  are  proposed  for  high-temperature  structural  applications  due  to  their  high-strength  to 
density  ratios  and  ability  to  withstand  high  temperatures.  High  temperature  capability  allows 
engines  to  operate  more  efficiently,  and  with  a  reduction  in  the  complex  cooling  systems  often 
required  for  metallic  alloys.  Performance  also  is  increased  by  weight  reductions  from  lighter 
materials.  Fiber  reinforcement  of  the  ceramic  matrix  allows  the  inherently  brittle  material  to  fail  in 
a  more  graceful  manner  due  to  crack  bridging  and  fiber  pull-out.  These  properties  make  ceramic 
matrix  composites  ideal  materials  for  inert  environments,  however,  carbon  fiber  reinforced  ceramic 
matrix  composites  are  prone  to  corrosion  in  oxidizing  environments. 

Many  studies  of  the  oxidation  of  carbon  fibers  and  pyrocarbon  interphases  in  ceramic 
matrix  composites  involve  thermogravimetric  analysis  (TGA)  in  which  weight  loss  of  the 
composite  is  monitored  over  time  [1, 2,  3,  4].  These  types  of  studies  can  provide  valuable 
information.  However,  the  oxidation  of  the  carbon  constituents  in  the  composite  for  unstressed 
and  stressed  states  is  very  different.  In  an  unstressed  state,  silica  forms  near  1200°C  and  cracks 
close  near  the  processing  temperature.  These  effects  protect  the  interior  of  the  composite  from  the 
outside  environment  [2, 4,  5].  However,  in  real  application  conditions,  wide  temperature  ranges 
and  stresses  will  prevent  cracks  from  sealing. 

The  oxidation  of  carbon  fibers  in  a  SiC  matrix  has  been  studied  experimentally  in  the 
stressed  state  [6].  Temperatures  ranged  from  350°C  to  1500°C  and  stresses  of  69  MPa  ( 10  ksi)  or 
172  MPa  (25  ksi)  were  applied.  Although  weight  could  not  be  monitored  as  in  TGA  experiments, 
the  times  and  strains  to  failure  of  the  composites  provided  an  indication  of  the  extent  of  oxidation. 
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Times  to  failure  suggested  two  regions  of  oxidation  kinetics  with  a  strong  dependence  on 
temperature.  In  the  temperature  range  of  350°C-750°C,  times  to  failure  dropped  sharply  as 
temperature  was  increased.  Samples  went  from  surviving  beyond  25  hr.  with  a  high  residual 
strength  to  failing  in  9 1  minutes  at  69  MPa  and  in  22  minutes  at  the  higher  stress  of  172  MPa  for 
temperatures  of  750°C.  In  the  higher  temperature  range  of  1000°C-  1500°C,  temperature 
differences  did  not  have  as  significant  an  effect  on  sample  life.  Sample  lives  remained  near  the 
same  level  as  those  at  750°C  with  some  deviations  in  times  to  failure  across  the  higher  temperature 
range. 

The  effects  of  oxidation  were  also  determined  from  analysis  of  polished  cross-sections 
and  fracture  surfaces  using  optical  and  scanning  electron  microscopes.  Oxidation  was  observed  to 
occur  preferentially  along  microcracks.  Polished  cross-sections  of  samples  tested  in  the  lower 
temperature  range  showed  oxidation  was  relatively  minimal  and  occurred  in  a  rough,  uneven 
pattern  near  the  edges,  but  oxidation  was  also  observed  deep  into  the  interior  of  the  section, 
suggesting  slow  oxidation  kinetics.  Sections  from  samples  tested  in  the  higher  temperature  range 
showed  oxidation  occurring  primarily  along  the  edges  in  a  uniform  manner  so  that  a  reaction  front 
could  be  seen.  This  gave  a  shrinking  core  type  of  effect  as  outer  fibers  were  consumed  as  the 
oxidation  process  moved  inward.  This  pattern  along  with  the  observation  that  fibers  deep  in  the 
interior  of  the  sample  remained  free  of  attack  suggests  the  oxidation  kinetics  are  much  faster  than 
in  the  lower  temperature  range.  This  change  in  the  rate  controlling  mechanism  from  the  lower 
temperature  regime  to  the  higher  temperature  regime  has  been  attributed  to  a  shift  from  the  kinetics 
being  reaction  rate  limited  to  being  diffusion  limited.  These  experimental  results  were  used  to 
provide  a  basis  for  developing  a  model  so  that  a  better  understanding  of  diffusion  controlled  and 
reaction  controlled  kinetics  could  be  gained. 

A  finite  difference  model  was  developed  which  simulates  tlie  diffusion  of  oxygen  into  a 
matrix  crack  bridged  by  an  array  of  carbon  fibers.  The  model  provides  a  visual,  qualitative  and 
quantitative  basis  for  understanding  the  oxidation  kinetics.  The  influence  of  important  variables: 
temperature  (T),  diffusion  coefficient  (D),  and  reaction  rate  constant  (K),  are  studied. 

By  obtaining  a  better  understanding  of  the  oxidation  kinetics  of  carbon  fibers  in  a  ceramic 
matrix,  the  best  approaches  for  protecting  against  oxidation  can  be  determined  depending  on  the 
application  conditions:  environment,  temperature,  stresses  and  required  component  lives.  One 
possible  approach  for  protecting  the  fibers  might  involve  applying  a  less  reactive  interphase  to  the 
fiber.  Another  may  be  to  include  an  oxidation  inhibitor  in  the  matrix  that  will  react  with  oxygen  to 
form  oxides  or  glasses  to  seal  cracks  or  to  consume  oxygen  to  provide  more  time  before  the  carbon 
is  attacked.  Still  another  approach  could  be  to  identify  external  seal  coatings  that  can  better  protect 
tlie  interior  of  the  material.  One  approach  may  be  preferred  over  another  due  to  exposure 
conditions,  cost  and/or  mission  requirements. 

THEORY 

The  related  theory  used  in  developing  the  finite  difference,  cracked  matrix  model  will  be 
briefly  discussed.  In  a  ceramic  matrix  composite,  oxygen  diffuses  down  a  pore  or  crack  and  reacts 
with  carbon  to  form  an  oxide  (either  carbon  monoxide  or  carbon  dioxide).  The  oxide  then  diffuses 
out  of  the  pore  to  the  atmosphere.  Non-Reaclive  Matrix 

See  Figure  1 . 

O2  +  C  — >  CO2  ^ 

02  +  2C  2CO 

Figure  1 .  Oxygen  diffusing  into  a  pore  and  reacting  with  carbon. 
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The  oxidation  kinetics  for  the  recession  of  carbon  can  be  described  by  the  linear- 
parabolic  rate  law: 

jr  +  x  =  t,  (1) 

kp  ^1  _  _  2/  \  i,  • 

where  x  (m)  is  the  recession  distance  of  the  carbon,  kp  is  the  parabolic  rate  constant  (m  /sec),  ki  is 
the  linear  rate  constant  (m/sec)  and  t  is  time  (sec). 

The  parabolic  and  linear  rate  constants  can  be  related  to  more  familiar  mass  transport 
constants,  the  diffusion  coefficient  and  the  reaction  rate  constant  respectively.  Only  the  relevant 
equations  will  be  discussed.  More  thorough  discussions  and  derivations  can  be  found  in  papers  by 
Eckel  et  al.  [7]  and  Glime  et  al.  [8]. 

When  the  linear  term  in  equation  (1)  dominates  and  carbon  dioxide  is  the  evolving  oxide, 
the  equation  for  the  recession  distance  is  expressed  as  [7]: 

x  =  ti/2  kp1/2  =  t1/2  (2D  Cox  /  pc )  such  that  kp  =  (2D  Cox  /  Pc )  (2) 

where  D  is  the  diffusion  coefficient  of  oxygen  in  the  pore  (m2  /  sec),  Cox  is  the  oxygen 
concentration  of  the  atmosphere  (mol  /  m'  )  and  pc  is  the  molar  density  of  carbon  (mol  /  m  ). 

In  the  case  where  carbon  monoxide  is  the  reaction  product.  Equation  (2)  for  the  recession 
distance  will  have  a  different  form.  However,  the  particular  oxide  that  forms  and  the  exact  value 
for  the  diffusion  coefficient  is  not  of  critical  importance  in  this  discussion  and  in  the  model.  The 
main  concern  at  this  point  in  the  development  of  the  model  is  to  identify  trends.  The  discussion 
does  not  consider  Knudsen  effects  in  which  molecular-wall  collisions  are  more  frequent  than  inter- 
molecular  collisions.  Knudsen  effects  are  not  a  factor  until  the  pore  becomes  as  small  as  0.1  pm 
[7].  In  the  material  studied,  typical  cracks  in  the  matrix  and  seal  coating  are  on  the  order  of  1  pm. 

The  second  limiting  case  occurs  when  kp  is  much  larger  than  k|.  The  recession  distance 
can  be  expressed  as  [8]: 

x  =  t  k,  =  t  K  Cox  /  pc  such  thatk,  =  K  Cox  /  pc  (3>- 

The  above  equation  takes  an  alternate  form  when  the  oxygen  concentration  is  substituted 
by  the  oxygen  partial  pressure, times  the  total  concentration  of  gas  molecules  in  the  atmosphere, 
CT  (mol/m3).  The  total  concentration  can  be  substituted  using  the  ideal  gas  law  so  that  Equation 
(3)  becomes  [7]: 

k,  =  ( 1/  pc  )(x  CT  )  (K)  =  (1/  pc)  ( X  P/RT)  K  (4) 

where  P  is  pressure  (Pa),  R  is  the  gas  constant  (J/mol  K),  and  T  is  the  absolute  temperature  (K). 
Substituting  values  for  kp  and  k| ,  Equation  (1)  becomes 

t  =  r  +  x  =  x2  +  _ x _  (5). 

kp  k|  (2D  Cox  /  pc )  (1/pcX  XP/RT)K 

The  diffusion  coefficient,  D,  can  be  obtained  from  Chapman-Enskog  kinetic  theory  [9]  and  the 
reaction  rate  constant.  K,  can  be  expressed  in  an  Arrhenius  form: 
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Dab  =  (5.9543  x  10'24)  T3/2(l  /Ma+\  /M_r)172  ,  K  =  ko  exp(-Q/RT)  (6)  &  (7) 

P  G~AB  ^ab 

where  DAB  (m2  /  sec)  is  the  diffusion  coefficient  of  gas  A  (oxygen)  diffusing  in  the  stationary  or 
stagnant  atmosphere  of  gas  B  (oxide),  T  is  temperature  (K),  MA  and  MB  are  the  molecular  weights 
of  gases  A  and  B  (kg/mol),  and  P  is  pressure  (Pa).  The  symbols  cAB  (m)  and  QAB  (dimensionless) 
are  Leonard  Jones  potentials  for  the  collision  diameter  and  collision  integral  respectively  .  Q  is  the 
activation  energy  (J/mol). 

From  the  above  equations  for  the  diffusion  coefficient  (6)  and  the  reaction  rate  constant 
(7),  it  is  seen  that  one  of  the  two  arguments  in  Equation  (5)  will  dominate  depending  on 
temperature.  The  diffusion  coefficient  depends  on  T3/2,  with  the  effect  being  even  less  when  the 
decreasing  molar  density  of  the  gas  is  considered  with  increasing  temperature.  The  reaction  rate 
constant  varies  exponentially  with  temperature  when  large  activation  energies  are  considered.  It 
can  be  seen  that  Equation  (5)  describes  a  two  step  process  occurring  in  series.  Each  step  has  the 
potential  to  limit  the  overall  process.  When  the  reaction  rate  constant  is  small  and  the  diffusion 
coefficient  is  large,  the  process  will  be  limited  by  the  reaction  rate,  i.e.  reaction  controlled. 

Oxygen  will  diffuse  into  the  pore  faster  than  it  can  react  so  that  the  pore  becomes  saturated  in 
oxygen.  However,  if  the  diffusion  coefficient  is  small  compared  to  the  reaction  rate  constant, 
carbon/oxygen  reactions  occur  as  soon  as  oxygen  is  supplied.  In  this  case  the  process  is  controlled 
by  the  slower  diffusion  step,  i.e.  diffusion  controlled. 

MODELING 

In  order  to  gain  a  better  understanding  of  the  diffusion  and  reaction  kinetics  observed  in 
the  experiments,  a  finite  difference  model  was  developed.  The  physical  setup  of  the  model  is  an 
array  of  carbon  fibers  bridging  a  matrix  crack,  as  illustrated  in  Figure  2.  Figure  3  shows  an 
example  of  an  experimental  observation  of  oxidation  occurring  along  fibers  bridging  a  crack, 
which  is  the  physical  situation  on  which  the  model  is  based. 


Figure  3.  Oxidation  of  carbon  fibers  that  bridge 
Figure  2.  Fibers  bridging  a  cracked  matrix.  matrix  cracks. 
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A  2-Dimensional  grid  pattern  was  laid  out  as  a  250x250  mesh  of  grids,  which  represents 
an  open  plane  within  the  bridged  space  parallel  to  the  two  crack  surfaces  of  the  matrix.  Grids 
within  the  mesh  are  designated  as  either  open  space  containing  the  diffusing  gases,  or  as  carbon 
containing  grids.  Carbon  fibers  are  represented  as  a  12x12  array  of  fibers  which  have  diameter  of 
10  grids.  Fibers  are  spaced  a  distance  of  10  grids  from  one  another  and  from  the  outer  edge.  One 
quarter  of  the  grid  pattern  is  shown  in  Figure  4. 


Midline 


Figure  4.  One  quarter  plot  of  the  matrix  pattern  (right)  orientated  as  shown  on  the  left. 

In  the  model,  the  outer  boundary  has  a  initial  oxygen  concentration  arbitrarily  set  at  1.0 
mol/cm3,  which  is  assumed  to  stay  constant  due  to  the  atmospheric  conditions.  Oxygen  diffuses 
into  the  matrix.  Concentrations  are  calculated  throughout  the  matrix  based  on  oxidation  kinetics  in 
which  both  reactivity  and  diffusion  are  both  dependent  on  temperature  and  one  may  be  limiting  the 
overall  process.  Plots  of  oxygen  concentration  as  a  function  of  grid  position  showing  gradients 
high  or  low  in  oxygen  concentration  will  give  clues  to  the  limiting  step  in  the  two  step  process. 

When  mass  transport  occurs  by  convection,  the  flux  is  given  by, 


Flux  =  J  =  K  (CA.i  -  C.O ) 


where  CA.i  (which  has  a  value  of  0  mol/cm3)  is  the  oxygen  concentration  at  the  reaction  front  or  of 
the  carbon  edge  (carbon  containing  grid)  and  C^  is  the  oxygen  concentration  in  the  bulk  gas 
(oxygen  containing  grid). 

When  mass  transport  occurs  by  diffusion,  the  flux  is  given  by, 

Flux  =  J  =  D  (C2  -  C]  )  /  Ax 


where  (C2  -  Cj  )  is  the  difference  in  oxygen  concentration  within  the  stagnant  gas  across  some 
distance,  Ax.  An  important  unitless  parameter,  the  Sherwood  number,  Sh,  is  obtained  from  the 
ratio  of  the  fluxes  per  the  driving  force  due  to  concentration  differences, 

Sh  =  F(convection)  /AC  /  F(diffusion)  /AC  =  K  x  /  D 

The  Sherwood  number  will  be  used  extensively  throughout  the  model.  For  low  Sh  values,  low 
temperatures,  the  kinetics  are  reaction  controlled,  and  for  high  Sh  values,  high  temperatures,  the 
kinetics  are  diffusion  controlled. 
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Detail  of  the  grid  pattern  showing  half  of  a  carbon  fiber  and  open  pore  space  is  illustrated 
in  Figure  5.  The  oxygen  concentration  of  a  grid  is  determined  by  sampling  from  its  four  neighbors. 
The  fluxes  between  the  central  grid  and  its  4  neighbors  are  added.  The  net  flux  into  and  out  of  the 
central  grid  is  zero.  The  flux  equation  is  then 
solved  for  the  concentration  of  the  central 
grid.  When  all  four  neighbors  are  within  the 
stagnant  gas  as  for  grid  1,  mass  transport 
occurs  by  diffusion  and  the  equation  for  the 
oxygen  concentration  of  this  grid  is 

C(i  j)  -  sum  of  oxygen  concentrations  (i,j-l)  [ 

4 

When  one  of  the  neighboring  grids  contains 
carbon,  as  for  grid  2  in  Figure  5,  mass 
transport  occurs  through  convection  between 
the  central  grid  and  the  carbon  containing 
grid.  For  the  other  three  grids  mass  transport 
occurs  by  diffusion.  The  equation  for  the 
oxygen  concentration  of  the  central  grid  becomes 

ca  j)  =  raij-iHCffKi)+aij+ni  +  sh*ai+ij) 

(3+Sh) 

When  two  neighbors  contain  carbon  and  two  contain  the  stagnant  gas,  as  in  grid  3,  the  equation 
becomes 

c(i  j)  =  rai.i-n+ai-i.m  +  SM£fij±ij±cfl±  i.m 

(2+2Sh) 

The  program  is  run  for  several  different  Sh  values:  0.0001, 0.001,  0.01, 0.1,  1,  10. 

Oxygen  concentrations  are  determined  throughout  the  bridged  section  based  on  the  boundary 
condition  and  the  equations  for  oxygen  concentration.  The  program  runs  until  steady  state  is 
reached.  Note  the  program  does  not  remove  carbon  by  reaction  at  this  point;  future  development 
of  the  program  will  remove  carbon  and  continuously  recalculate  oxygen  concentrations  so  that 
oxidation  patterns  relating  to  diffusion  controlled  and  reaction  controlled  kinetics  are  obtained. 

The  program  may  be  thought  of  at  this  point  as  having  a  continuous  supply  of  oxygen  and  carbon, 
where  oxygen  concentrations  throughout  the  section  are  calculated  until  steady  state  is  reached.  It 
should  be  pointed  out  that  since  oxidation  really  occurs  in  3-dimensions,  the  2-  dimensional  model 
does  not  consider  or  allow  the  analysis  of  certain  factors.  These  include  3rd  dimensional  effects 
such  as  fiber  recession  and  the  supply  of  oxygen  through  open  channels  where  the  carbon  fiber  has 
been  consumed. 

The  obtained  data  files  for  a  given  Sherwood  number  were  used  to  prepare  plots  of  the 
oxygen  concentrations  across  the  interior  of  the  fiber  bridged  matrix.  Gradients  in  oxygen 
concentration  across  the  bridged  section  show  regions  where  oxygen  supply  is  low  or  high  and 
give  indications  about  the  reactivity.  Figure  6  shows  the  oxygen  concentrations  for  one-quarter  of 
the  whole  section  for  a  relatively  high  Sherwood  number,  Sh=0.1.  The  3-D  effect  of  the  plot  is  due 
to  plotting  oxygen  concentrations  across  the  2-D  surface.  It  is  seen  that  oxygen  concentrations 
drop  sharply  upon  diffusion  into  the  matrix.  This  suggests  that  the  carbon  is  very  reactive  and 
reacts  with  oxygen  as  soon  as  it  is  supplied.  The  overall  oxidation  process  is  limited  by  the  slow 
supply  of  oxygen  to  the  reaction  process  at  the  carbon/oxygen  interface,  i.e.  diffusion  controlled 
kinetics.  The  case  for  a  relatively  low  Sherwood  number,  Sh=0.001,  is  shown  in  Figure  7.  The 
plot  shows  that  the  reactions  are  slow  so  that  oxygen  is  able  to  diffuse  into  the  interior  in  high 
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concentrations  and  saturate  the  interior.  In  this  case,  the  two  step  process  of  oxygen  diffusion  and 
carbon/oxygen  reactions  is  limited  by  the  slower  reaction  process,  reaction  controlled  kinetics. 


Figure  6  (left)  Sh=0.1.  Figure  7  (right)  Sh=0.001.  One -quarter  plots  of  a  250x250  matrix  with  a  12x12  fiber 
array.  Interated  to  equilibrium. 

In  the  two  previous  figures,  the  curves  at  the  edges  of  the  plots  represent  the  oxygen 
concentration  through  one -half  of  the  thickness  of  the  section.  Full  plots  of  the  oxygen 
concentrations  across  the  midline  are  shown  in  Figure  8  for  several  Sherwood  numbers.  High  Sh 
values  illustrate  how  the  process  is  diffusion  controlled.  Oxygen  is  quickly  consumed  before  it  can 
enter  into  the  interior  of  the  matrix.  Low  Sh  values  illustrate  how  the  kinetics  are  reaction 
controlled.  Carbon/oxygen  reactions  are  slow  so  that  there  is  a  high  concentration  of  oxygen 
throughout  the  section. 


Sh  =  Ax  K 
D 

Reaction  Controlled 


Sh  =  1 

and  10  | 

(lower) 

Diffusion  Controlled 


0  50  100  150  200  250 

Position  Across  Matrix  (0-250) 

Figure  8.  Oxygen  concentrations  across  the  center  of  the  surface  for  several  Sherwood  (Sh)  values. 
SUMMARY 

Experimental  studies  on  C/SiC  composite  material  exposed  to  stressed  oxidation 
conditions  suggested  two  types  of  oxidation  kinetics.  Lives  of  the  samples  were  very  temperature 
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dependent  in  the  range  of  350°C-750°C,  while  there  was  less  of  a  temperature  dependence  in  the 
range  of  1000°C-1500°C.  Microstructural  analysis  revealed  a  different  type  of  oxidation  pattern 
for  each  temperature  range.  The  related  theories  for  carbon  recession,  reaction  kinetics,  diffusion 
kinetics,  flux  due  to  mass  transport,  and  Sherwood  numbers  were  reviewed.  Temperature  effects 
on  the  reaction  rate  constant  and  the  diffusion  coefficient  were  also  discussed.  Based  on 
experimental  observation,  and  theory,  a  finite  difference  model  was  presented. 

CONCLUSIONS 

The  oxidation  of  carbon  fibers  in  a  ceramic  matrix  is  a  two  step  process  with  a  varying 
temperature  dependence.  The  overall  process  will  be  limited  by  the  slower  step.  The  kinetics  can 
be  analyzed  by  studying  the  diffusion  of  oxygen  and  determining  its  concentration  throughout  the 
section.  The  finite  difference  model  is  in  good  agreement  with  experimental  results  in  suggesting 
two  types  of  kinetics.  The  model  supports  the  conclusion  that  the  kinetics  are  reaction  controlled 
at  low  temperatures,  and  diffusion  controlled  at  high  temperatures. 

The  next  steps  in  the  model  development  include  the  removal  of  carbon  as  it  reacts  with 
oxygen.  In  the  case  of  diffusion  controlled  kinetics,  i.e.  high  Sherwood  number,  oxygen  will 
oxidize  away  the  outer  fibers  first,  then  continue  to  move  inward  so  that  a  shrinking  core  effect  is 
seen.  For  reaction  controlled  kinetics,  i.e.  low  Sherwood  number,  there  will  be  a  generalized 
oxidation  of  fibers  throughout  the  section.  Another  developmental  step  for  the  model  will  include 
putting  an  element  around  the  fibers  to  represent  the  interphase.  Different  reactivities  of  the  fiber 
and  interphase  can  then  be  considered. 
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